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^ Abstract 

The current constraints from lattice QCD on the existence of the H-dibaryon are discussed. With 
only two significant lattice QCD calculations of the H-dibaryon binding energy at approximately 
the same lattice spacing, the forms of the chiral and continuum extrapolations to the physical point 
are not determined. In this brief report, we consider the constraints on the H-dibaryon imposed 
by two simple chiral extrapolations. In both instances, the extrapolation to the physical pion 
mass allows for a bound H-dibaryon or a near-threshold scattering state. Further lattice QCD 



calculations are required to clarify this situation. 
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The possibility of a bound H-dibaryon (a six-quark state uuddss) was first proposed by 
Jaffe [lj in 1977 using an MIT bag-model |2j calculation which predicted that such a state 
would have a mass 81 MeV below the AA threshold. Since then, experiments [3HH] and other 
theoretical calculations [7] |8], including a few pioneering lattice QCD calculations [HHI1], 
have been seeking evidence of its existence. In 1997, Bashinsky and Jaffe [T5] reviewed 28 
model predictions on the mass of the H-dibaryon, which vary from 1.1 to 2.9 GeV. (See 
Fig. 4 in Ref. [15] for a histogram of model predictions.) Two years later, Sakai et al. [7j 
comprehensively reviewed various updated experimental and theoretical results (see Fig. 1 
in Ref. [7j), extending the upper limit of model predictions to around 3.1 GeV. Model 
predictions for the H-dibaryon mass have not changed much since then. Recently, the 
NPLQCD and HALQCD collaborations reported results showing that the H-dibaryon is 
bound for a range of light-quark masses larger than those found in nature [T6"| [T7] . These 
calculations are important for a number of reasons. First, they show that lattice QCD is 
now capable of calculating the energy of simple nuclei, with the H-dibaryon being an exotic 
example. Second, they provide evidence that a bound H-dibaryon may exist for some values 
of parameters entering the QCD Lagrangian. However, it is important to determine if this 
system is, in fact, bound at the physical values of the light-quark masses and with the 
inclusion of electroweak interactions. Experimental evidence currently suggests that such 
a bound state does not exist between 2.136 and 2.231 GeV [3] (indicating either that the 
binding energy could be larger than 95 MeV or there is no bound H state at all), but that 
a near-threshold resonance may exist in the scattering channel with the quantum numbers 
of the H-dibaryon [4.J. In this note we investigate the current constraints on the binding 
energy of the H-dibaryon at the physical value of the pion mass, in the isospin limit and in 
the absence of electroweak interactions, by extrapolating the available lattice QCD results. 

Early quenched lattice QCD calculations [91414"] drew inconsistent conclusions regarding 
the binding of the H-dibaryon. Table [I] shows a summary of parameters used in those works 
and their conclusions. Generally the calculations involved less than 100 measurements, and 
the analyses lacked estimates of systematic uncertainties, such as 0(b) discretization errors 
for Wilson fermions, finite- volume effects, and quenching (could be 30% or larger depending 
on the physical quantity). 

TABLE I: Summary of H-dibaryon lattice QCD calculations. The pion masses used in Refs. [X 11 — 
[L3] are estimated from the bare quark-mass parameters listed in their proceedings. For the works 
that do not perform volume extrapolation to the V — > oo limit, the mass difference between the 
H-dibaryon and 2-A state from the largest volume is used. 



Ref. 


n f 


Sf 


a' 1 (GeV) 


m n (GeV) 


L (fm) 


V -> oo 


m H - 2m A (MeV) 


Mackenzie et al. [9] 





Wilson 


0.9 


0.60-0.90 


1.3 


no 


> 


Iwasaki et al. [10] 





Wilson 


1.81 


0.70-1.90 


1.74 


no 


[-676,-381] 


Pochinsky et al. [TT] 





Wilson 


1.79 


> 0.5 


2.08, 3.12 


yes 


> 


Wetzorke et al. [TJ [T3] 





clover 


1.11 


> 0.5 


1.4-4.27 


no 


> 


Luo et al. [H] 





aClover* 


1.48 


0.70-0.85 


3.2,6.4 


no 


-64(59) 


NPLQCD (this work) 


2+1 


aClover* 


5.6 


0.39 


2.0, 2.5, 3.0, 4.0 


yes 


-13.2(1.8)(4.0) 


NPLQCD (this work) 


2+1 


aClover* 


5.6 


0.23 


4.0 


no 


0.6(8.9)(10.3) 


HALQCD [T7] 


3 


clover 


1.63 


0.67-1.01 


2.0, 3.0, 4.0 


yes 


[-40, -30] 



* The abbreviation "aClover" stands for the anisotropic clover fermion action. 
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Unlike the single-hadron spectroscopy of a typical lattice QCD calculation, high statistics 
are essential to manage the uncertainties in the extraction of bound-state information. To 
demonstrate this, we have broken down our measurements into smaller subsets to see how 
the uncertainty evolves [IS]. Figure 28 in Ref. [TS] shows NPLQCD's extracted values of 
(fccot5) _1 for AA scattering (5: phase shift) versus the extracted value of |k| 2 (k: center- 
of-mass momentum) with total number of measurements varying from 10 4 to 4 x 10 5 in a 
lattice with spatial extent of 2.5 fm. With a small number of measurements, one could easily 
extract a value of k cot 5 that appears repulsive at this particular volume without a careful 
estimation of the systematic uncertainty. Even with 10 5 total measurements, we are not 
yet able to conclusively determine whether the interaction is attractive or repulsive. It is 
only when the total number of measurements approaches 4 x 10 5 that the interaction can 
be determined to be attractive, and even then with only about 2a significance. For the 
remainder of this brief report, we will concentrate on the two recent dynamical studies with 
higher statistics by the NPLQCD and HALQCD collaborations. 

The details of the two lattice QCD calculations that provide statistically significant ev- 
idence for a bound H-dibaryon can be found in the recent works of NPLQCD [TB] and 
HALQCD [17] . The NPLQCD result is determined from nf = 2 + 1 calculations in four 
lattice volumes (with spatial extents of L ~ 2.0, 2.5, 3.0 and 4.0 fm), each at a single spatial 
lattice-spacing of b ~ 0.123 fm and a pion mass of ~ 390 MeV. Ref. [16] reported a bind- 
ing energy of Bh = 16.6 ± 2.1 ± 4.6 MeV at that pion mass; with more statistics on this 
ensemble, the updated binding energy is now 13.2 ± 1.8 ± 4.0 MeV. The HALQCD collab- 
oration performed calculations in three lattice volumes (with spatial extents of L ~ 2.0, 3.0 
and 4.0 fm) at a lattice spacing of b ~ 0.121 fm and in the limit of SU(3) flavor symmetry 
(rif = 3) at three different quark masses giving m n ~ 673,837 and 1015 MeV. In order to 
extrapolate in the light-quark masses, the binding energy of Bh = 37 .4 ± 4.4 ± 7.3 MeV 
obtained at m n ~ 837 MeV is used because this pion mass corresponds to a strange-quark 
mass that is closest to that of nature 1 (and to that of the NPLQCD calculations) 2 . 

NPLQCD and HALQCD employed different clover discretizations of the light quarks (and 
different gauge-actions), providing results that are (9(6)-improved and therefore both sets 
of calculations have lattice-spacing artifacts that scale as 0(b 2 ). Given the precision with 
which the single-hadron energy-momentum relation is satisfied [IB] , the contributions from 



At leading order in chiral perturbation theory = 2Am and m 2 K — A(m s + m) (Gell-Mann-Oakcs- 
Renner) in the limit of isospin symmetry, where m is the isospin-averaged light quark mass and A is related 
to the quark condensate. Inserting the NPLQCD values of the pion and kaon masses into this relation 
would lead to using the H-dibaryon binding energy calculated at ~ 673 MeV by HALQCD. However, 
lattice QCD calculations (in particular, calculations on additional ensembles that have the same strange 
quark mass as those used in the NPLQCD calculations, as shown in Table VI of Ref. [IH]), which include 
the full light-quark mass dependence, dictate that the binding energy calculated at m n ~ 837 MeV is the 
appropriate value to use for the SU(3) symmetric point. Given the size of the uncertainties of both lattice 
QCD calculations, the difference in the results obtained with the ~ 837 MeV and m n ~ 673 MeV are 
not statistically significant. 

In this channel, the ground-state energy-levels determined in the lattice QCD calculations of 
HALQCD [T7] are expected to be exponentially close to the actual bound-state energy of the H-dibaryon 
and do not suffer from the uncontrolled approximations that are present [5D] in phase-shifts calculated 
via the energy-dependent and sink-dependent potentials presented by HALQCD. 
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Lorentz-symmetry breaking operators that appear at order 0(b 2 ) are expected to be highly 
suppressed. Naive scaling arguments, as well as the cancellations that occur in forming 
energy differences, suggest that lattice spacing artifacts are suppressed, as compared, for 
instance, to the leading quark-mass effects. However, definitive statements about the lattice 
spacing dependence will require calculations at smaller lattice spacings. In this context, 
we assume that lattice spacing artifacts are small (and therefore are not driving the H- 
dibaryon binding), and focus on the light-quark mass dependence. Alternately, one can 
view the present study as performing the chiral extrapolation of the (approximately) fixed 
lattice-spacing results. 

The form of the chiral extrapolation of the H-dibaryon binding energy is unknown, and 
there are a number of features that make it difficult to determine. In the limit of exact 
SU(3) flavor symmetry, the non-interacting ground state of the I = 0, J = and s = —2 
(the quantum numbers of AA) system is multiply degenerate, comprised of the hadronic 
states AA, EE, and AS. The interactions among these states produce the H-dibaryon 
(when it is bound) as the ground-state, with the other states being of higher energy. Even 
in the presence of SU(3) breaking, as is the case in the NPLQCD calculations, the spectrum 
is qualitatively similar with the non-interacting states being nearly degenerate. If the H- 
dibaryon is tightly bound, then a chiral expansion 3 of the form of that for single hadrons will 
result. On the other hand, if the H-dibaryon is loosely bound, the binding energy results 
from a cancellation between short-distance and long-distance contributions, as is the case 
with the deuteron where the extrapolation has significant structure [2HH2E]- However, in 
contrast with the NN system, the long-distance contribution to AA interactions results from 
two-pion exchange. Consequently, the H-dibaryon is expected to be less fine tuned than the 
deuteron, leading to the expectation that the chiral extrapolation may resemble the form 
B H (m n ) = B + d x m\ + C(m 3 ). 

The quark mass dependence of the lowest-lying octet baryon masses, such as the N, A, 
E and S, is an ongoing topic of discussion, and motivates a second extrapolation form. The 
form of the light-quark mass dependence of the baryon masses produced by the lowest orders 
of baryon chiral perturbation theory does not naturally reproduce the currently available 
results of lattice QCD calculations 4 . This is currently interpreted as an indication that 
the chiral expansion is not converging at these heavier pion masses and the possibility 
there is a new scale associated with the strong interactions [30]. At the least, there is 
correlation between higher orders in the expansion beyond the expectations from naive 
dimensional analysis. Walker-Loud [31] performed a comprehensive analysis of all of the 
lattice QCD calculations of the nucleon mass and found that the results are consistent 
with linear dependence on the dimensionless variable m n / fi°\ where /i is the pion decay 
constant in the chiral limit, with an extrapolation that is consistent with the experimental 
value. This is in contrast with the expectations of an expansion about the chiral limit, which 
has the form Mjj{m^) = M + a 2 fn^ + 0(m^), where M and are parameters that 
must be determined from the lattice QCD calculations. This motivates us to consider a 
non-analytic extrapolation of the form Bnim^) = Bq + cim n + 0(m 2 ). It is possible that 
the true form lies somewhere between the linear and the quadratic forms, with cancellations 



3 Low-energy effective field theories have been constructed to describe baryon systems carrying strangeness; 
see Refs. [SH24]. 

4 SU(2) heavy-baryon chiral perturbation theory does reproduce the lattice results for the nucleon mass 
(with x 2 /dof ~ 1), but with large cancellations between different orders in the expansion [29 . 
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occurring between analytic and non-analytic contributions. Without any better guidance as 
to the form of the chiral extrapolation, we will consider the results from these two forms of 
extrapolation with relatively heavy pions to provide nothing more than an estimate of the 
H-dibaryon binding energy at the physical light-quark masses. Extrapolations with more 
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FIG. 1: The results of lattice QCD from the NPLQCD collaboration |16| and this work (n/ = 
2 + 1) (red squares), and the HALQCD collaboration [T7] (rt/ = 3) (blue triangles). The filled 
symbols are used in the extrapolations, while the open squares (NPLQCD's 230-MeV data) are 
not. Left panel: The darker (lighter) shaded region corresponds to an extrapolation of the LQCD 
calculations that is quadratic in the pion mass, of the form Bu^^) = Bq + d\ where 
the parameters are determined by the central values and statistical uncertainties (statistical and 
systematic uncertainties combined in quadrature). The vertical dashed (green) line corresponds to 
the physical pion mass. Right panel: Same as the left panel but with the extrapolation Bnim^) = 
B + ci m^. 

complicated behaviors are allowed but cannot be constrained by the current lattice QCD 
calculations, and are not discussed further. The chiral extrapolation of the H-dibaryon 
binding energy using the form Bnifn^) = Bq + d± results in the shaded region shown in 
Fig. [I] (left panel). The H-dibaryon binding energy at the physical value of the pion mass, 
neglecting isospin-violation and electromagnetic interactions, is found to be 

£^ adratic = 7.4 ±2.1 ±5.8 MeV , (1) 

as indicated by the intercept of the shaded region with the (green) dashed line in Fig. [I] (left 
panel). The first uncertainty results from an extrapolation using the statistical uncertainties 
of both lattice QCD calculations, while the second uncertainty results from the systematic 
uncertainties. The quadratic extrapolation suggests that the H-dibaryon is bound at the 
physical value of the pion mass. However, the H-dibaryon is unbound at the 2a level, and 
a near threshold scattering state remains allowed by the current lattice QCD calculations. 
Further, at the 2a level, the extrapolation is also consistent with the binding energy being 
independent of m v . 

Using the form Buim^) — Bq ± c\ to chirally extrapolate the lattice QCD calculations 
produces the results shown in Fig. [I] (right panel). The H-dibaryon binding energy at the 
physical value of the pion mass is found to be 

^lmear = _ Q 2 ± 3 3 ± 7 3 Me y ^ ^ 

as indicated by the intercept of the shaded region with the (green) dashed line in Fig. [I] (right 
panel). With the precision of the current lattice QCD results, the linear chiral extrapolation 
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does not discriminate between a bound or unbound H-dibaryon at the physical pion mass. 
However, it suggests that there is a state in the I = 0, J = and s = —2 channel that is 
either just bound or just unbound. 

Currently, both NPLQCD and HALQCD are improving their calculations to include 
lighter pion masses. NPLQCD is working on a calculation using annj = 2+1 ~ 230 MeV 
ensemble, which at sufficiently high statistics would improve our ability to extrapolate the 
H-dibaryon binding energy to the physical limit and determine whether it remains bound. 
To date, a single lattice volume with spatial extent 4 fm (giving m n L around 4.7) has 
been completed with around 9 x 10 4 measurements performed, and a binding energy of 
Bh = —0.6 ± 8.9 ± 10.3 MeV is obtained. This point is included for illustrative purposes in 
Fig. [TJ shown as an open square in both figures. However, an infinite volume extrapolation is 
currently unavailable at this pion mass, and we have not estimated the associated systematic 
uncertainty. Unfortunately, at this stage, this additional information does not differentiate 
between the extrapolation forms, and further improvements are needed before this light 
pion-mass ensemble can provide significant input to the extrapolation. 

We conclude that the current lattice QCD calculations of NPLQCD and HALQCD chi- 
rally extrapolated with plausible forms for the light-quark mass dependence, are not suffi- 
ciently precise at light enough quark masses to determine whether QCD predicts a bound 
H-dibaryon. However, the current lattice QCD calculations hint that the H-dibaryon be- 
comes less bound as the light-quarks become lighter, and the extrapolations suggest that 
there is a near-threshold state in the spectrum. Calculations at smaller lattice spacings 
are also required to ensure that lattice artifacts are not driving the observed H-dibaryon 
binding, and to allow extrapolation to the continuum limit. 

This result should be viewed in the context of the present experimental constraints, which, 
at face-value, effectively eliminate the possibility of a bound H-dibaryon [3] 5 . However the 
suggestion of structure in the scattering amplitude near threshold [4J would not be incon- 
sistent with the current lattice QCD results. In order to refine the lattice-QCD prediction, 
precision calculations at lighter quark masses are required. We are currently increasing the 
statistics of the calculations on the 230-MeV ensemble; this would significantly improve our 
ability to predict whether the H-dibaryon is bound or unbound at the physical light-quark 
masses. A multi-lattice-spacing, multi-volume lattice QCD study at the physical pion mass 
would be unambiguous, up to systematic uncertainties that arise from working in the isospin 
limit in the absence of electromagnetism. Ultimately, these effects will also be included. 

We would like to thank R. Jaffe for prompting us to consider the chiral-extrapolation, 
and D. Hertzog for emphasizing the importance of performing the extrapolation. 
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